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ABSTRACT 

We study the effect of large-scale magnetic fields on the non-axisymmetric Rossby 
wave instability (RWI) in accretion discs. The instability develops around a density 
bump, which is likely present in the transition region between the active zone and dead 
zone of protoplanetary discs. Previous works suggest that the vortices resulting from 
the RWI may facilitate planetesimal formation and angular momentum transport. We 
consider discs threaded by a large-scale poloidal magnetic field, with a radial field 
component at the disc surface. Such field configurations may lead to the production 
of magnetic winds or jets. In general, the magnetic field can affect the RWI even when 
it is sub-thermal (plasma /? 10). For infinitely thin discs, the instability can be 
enhanced by about 10 percent. For discs with finite thickness, with a radial gradient 
of the magnetic field strength, the RWI growth rate can increase significantly (by a 
factor of ^ 2) as the field approaches equipartition (/3 ^ 1). Our result suggests that 
the RWI can continue to operate in discs that produce magnetic winds. 

Key words: accretion, accretion discs - hydrodynamics - waves - planet formation - 
magnetic field 



1 INTRODUCTION 

Radial density structures can be common in astrophysical 
discs. For instance, in proto-planetary discs, different de- 
grees of coupling of magnetic field with disc gas could lead 
to a transition region where the disc material piles up (e.g., 
Gammie 1996; Matsumura & Pudritz 2003; Terquem 2008). 
Proto-planets can also induce dips/gaps in the disc through 
tidal interaction (Lin & Papaloizou 1986; Ward 1997). In 
accretion discs around black holes or neutron stars, sharp 
density contrasts can be created near inner disc edge in- 
side which the disc gas plunges inwards (e.g. Lovelace et 
al. 2009). 

Such density structures lead to shear flows that are 
susceptible to Kelvin-Helmholtz instability (Papaloizou & 
Pringle 1985; Papaloizou & Lin 1989; Lithwick 2009). 
A specific realization of this instability is the so-called 
Rossby Wave Instability (RWI) (Lovelace et al. 1999; Li et 
al. 2000,2001). For two-dimensional (vertically integrated) 
barotropic discs, the RWI relies on the existence of an 
extremum in the background fluid vortensity (also called 
potential vorticity; see Drazin & Reid 2004; Narayan et 
al. 1987), defined by C = k^/{2Q.Y.), where E is the disc 
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surface density, r2 the disc rotation rate and = AQ.^ + 
2rQdQ,/dr is the square of the radial epicyclic frequency. 
Since Rossby waves propagate along the gradient of vorten- 
sity, the instability can be understood as arising from the 
interaction between two Rossby waves propagating on each 
side of the vortensity extremum (e.g.. Tagger 2001; Lai & 
Tsang 2009). 

The nonlinear development of the RWI in discs leads to 
the production of vortices, which can accelerate the growth 
of meter-sized solids by concentrating dust grains and thus 
enhance the formation of planetesimals (Barge & Sommeria 
f 995; Johansen et al. 2004; Varniere & Tagger 2006; Heng & 
Kenyon 20f0). Such vortices may also induce sizable angular 
momentum transfer in the radial direction (Li et al. 200f). 
The vortices formed at the edge of planetary gaps can affect 
the rate of planet migration (de Val-Borro et al. 2007; Yu et 
al 20f0; Lin & Papaloizou 20ifa,b). 

While most studies of the RWI deal with two- 
dimensional (height-integrated) discs, Meheut et 
al. (2010, 20f2) carried out the global simulations of 
the RWI in 3D discs [see also Lyra & Mac Low (2012) for 
MHD simulations of 3D discs with no vertical stratifica- 
tion]. They showed that the Rossby vortices have significant 
vertical structure, which affects the long-term evolution 
of the vortices. Meheut, Yu & Lai (20i2) and Lin (20i2) 
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studied the linear RWI in 3D vertically stratified discs, 
showing that while the RWI is basically a 2D instability, 
appreciable vertical velocities can be induced in the disc. 

Magnetic fields are present in most accretion discs and 
can greatly change the dynamics of discs. It is generally 
thought that turbulence induced by the magneto-rotational 
instability (MRI; Balbus & Hawley 1998) is responsible for 
angular moementum transport in ionized discs. Recent simu- 
lations including ambipolar diflusion, however, suggest that 
the MRI associated with small-scale magnetic field may not 
be adequate to account for the angular momentum trans- 
port in protoplanetary discs, even in the "active zone" (Bai 
& Stone 2011). On the other hand, large-scale poloidal mag- 
netic fields threading the accretion disc can generate out- 
flows/jets and play an important role in the angular momen- 
tum transport (e.g., Blandford & Payne 1982). Large-scale 
magnetic fields are naturally present in the early phase of 
star formation. They may be advected inward with the ac- 
cretion fiow (e.g., Lubow et al. 1994), building up significant 
strength in the disc. The large-scale magnetic field may also 
counteract the disk self-gravity effect, thereby infiuencing 
disc fragmentation (Lizano et al. 2010). 

Recently, Yu & Li (2009) showed from linear analysis 
that a weak (sub-thermal) toroidal magnetic field in the disc 
can completely surpress the RWI. This surpression results 
from the modification by the magnetic field to the wave ab- 
sorption at the corotation resonance (see Fu & Lai 2011). In 
this paper, we study the effects of large-scale poloidal mag- 
netic fields on the RWI. As noted above, such a large-scale 
field is essential for producing bipolar outflows and may be 
necessary for angular momentum transport in the disc. Since 
the growth of the RWI is primarily due to in-disc motion 
(Meheut et al. 2012), we will focus on 2D disc dynamics in 
this paper. However, we will consider the effect of flnite disc 
thickness on the RWI growth rate while neglecting vertical 
stratification. 

Note that in our models (see Section 2), the large-scale 
poloidal magnetic field provides a coupling between the dy- 
namics of the disc and the magnetosphere (assumed to be 
current-free). We consider perturbations that have no ver- 
tical structure inside the disc (i.e., the vertical wavenum- 
ber kz = 0). Thus, we do not include MRI, which gener- 
ally involves perturbations with finite k^. Our basic disc- 
magnetosphere setup is self-consistent. Similar setups have 
been considered by various authors in different contexts 
(e.g.. Spruit et al. 1995; Tagger & Pallet 1999; Lizano et 
al. 2010). 

Our paper is organized as follows. In §2, the basic disc 
equations with large-scale magnetic fields are derived. In §3, 
we present the results of RWI of thin magnetized discs. In 
§4, we consider the effects of flnite disc thickness, and we 
conclude in 55. 



2 BASIC EQUATIONS 

We consider a thin conducting disc threaded by a magnetic 
field. The height-integrated mass continuity equation and 



momentum equation read 
r)y 

^ + Vx-(Su) = 0, (1) 

5-4^-^+4^^^^]-+^' (2) 

where V_l is 2D operator (acting on the disc plane), E, P 
and u are the surface density, height-integrated pressure and 
height-averaged velocity, respectively, g = —g{r)f (with g = 
GM/r'^) is the gravitational acceleration, [B]l = B(2: = 
H) — 'B{z = —H) (with H the half-thickness of the disc), and 
we have used [-B'^]l ~ 0. In this section we are considering 
infinitely thin discs, i.e., in the limit H/r — >■ 0. As a result, 
we neglect the internal magnetic force acting on the disc 
in Equation ((2)1. The effect of finite disc thickness will be 
considered in Section 4. 

For an equilibrium disc, the unperturbed velocity is 
u = (0, rn,0) (in cylindrical coordinates), with the angu- 
lar velocity determined by 



Q.r = 



(3) 



where = Br{z = H) — —B~. The unperturbed disc has 
-Sj — 0. Note that Br is nonzero only outside the discs, 
which has different signs just above and below the disc. In- 
side the disc, Br is zero and the differential rotation of the 
disc will not lead to generation of inside the disc. 

Note that the disc dynamics are coupled with the 
large scale poloidal magnetic field outside the disc, i.e., the 
disc magnetosphere. Now consider small-amplitude pertur- 
bations of the disc. Assuming that all perturbed quantities 
are proportional to exp(im</> — iu)t), where m = 1, 2, ■ ■ ■ is 
the azimuthal wave number and uj is the complex frequency, 
then the linearized fluid perturbation equations read 

= -Vx • (E Sm) , (4) 

d B 
—iuj Sur — 2n Sudi = — Sh + -— ^ SBr 
or 27rE 

where iIi — lj~ mO, is the Doppler-shifted frequency, k is the 
radial epicyclic frequency, and Sh = SP/Yl is the enthalpy 
perturbation (we assume barotropic discs). The magnetic 
field induction equation in the disc reads 



- iCj 5B^ = -Vx ■ {Bz 5u) , 
Combining Eqs. Q and 0, we have 



-6 



dr 



(7) 



(8) 



where = 5ur/{—iuj) is the Lagrangian displacement. 
Equation ((8]) can also be derived from {d/dt){Bz/'E) — 0. 
In terms of and Sh, the magnetic field perturbation is 



SB, = DiCr + DiSh, 



with 



c 



(9) 
(10) 



where Cs is the sound speed. 



Rossby Wave Instability in Magnetic Discs 3 



To determine SB:^ and SB'^ , we assume that the mag- 
netic field outside the disc is a potential field (e.g., Spruit et 
al. 1995; Tagger & Pallet 1999). This amounts to assuming 
that the Alfven speed above and below the disk is suffi- 
ciently high that currents are dissipated quickly (on the disc 
dynamical timescale). Define the magnetic potential S^m 
outside the disc via 



6B = -sign(2)V(5$j, 



(11) 



Then 5^m satisfies the Poisson equation (Tagger & Pallet 
1999) 



~2 5B,S{z) 



(12) 



where 5{z) is the Dirac delta function. The integral solution 
of (HU is 

S<fM{r) = I 5B4r') [| 6r/2(«)] dr' , (13) 
where a — r' /r, and the Laplace coefficient is defined by 



(14) 



TT Jq (a^ + Eq -I- 1 — 2q cos ( 

with eo the softening parameter (of order the dimensionless 
disk thickness H/r). The perturbed radial magnetic field at 
the upper disc surface is 

d 



sb: 



dr 



lsB^{r')[ 



d&r/2(Q) 
da 



dr', (15) 



and the azimuthal field is 5B^ — —{im/r)5$M ■ 

In our numerical calculations, we use the variable 
and Sh. The corresponding perturbation equations are 



d^r 
dr 



dr 
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Sh + 
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r^uj'-^ 2itT, 



-2 2^ B+ \ 2mn 
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27rE / ru! 

B^ dS^M 2mn B, „ 
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27rE dr 



rCo 27rE 



(16) 



(17) 



3 RESULTS FOR THIN DISCS 

We consider a density bump in the disc of the form 



S(r) 
Eo 



1 -f- (A - 1) exp 



(r- 



ro) 



2A2 



(18) 



We adopt A — 1.2 and A = 0.05ro throughout this paper. 
We further assume Cs — O.lroSlo ~ constant. 

In this section, Bz is taken as a constant independent 
of r (we will consider varying Bz in Sections 4-5) . The mag- 
nitude of Bz is specified by the dimensionless ratio 

^ ~ /v^ o2 vr72' (1*^) 
[T.oil^roy'^ 

where the subscript "0" implies that the quantities evaluated 
at r = ro. The corresponding plasma /3 parameter in the disc 
at r = ro is 



~Br 



An Ho '2 
Bz , 

ro 



(20) 
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Figure 1. Vortensity profiles for magnetized discs with Bz = 
0.78, corresponding to the plasma /9 ~ 2. The solid line is for the 
case of = and the dashed line for B^ = Bz ■ 



where we have used H = Cs/0.k — Cs/^l. We solve for the 
equilibrium rotation profile using Eq. (|3]). In Fig. [T] we show 
two examples of the vortensity profiles for the B^ = and 
B^ = Bz cases. Note that when B^ = 0, the vortensity pro- 
file is the same as a nonmagnetic disc. If B^ ^ 0, when Bz 
increases, the epicyclic frequency approaches zero gradually 
at the minimum, as does the vortensity. 

Since the eigenfrequency uj is in general complex, equa- 
tions (I16p - (|17p are a pair of first-order differential equa- 
tions with complex coefficients which are functions of r. We 
solve these equations using the relaxation method (Press et 
al. 1992), replacing the ODEs by finite-difference equations 
on a mesh of points covering the domain of interest (typ- 
ically 0.4 < r/r-o < 1.6). According to Eqs. dSJ, dTS]), and 
P5|l . both the magnetic potential and its derivative can be 
expressed in terms of linear combination of and Sh. For 
numerical convenience, we calculate the terms in the square 
bracket in Eqs. (|13p and p5|l and store them for later use. 
Note that these terms are computed only once and can be 
used repetitively. The wave equations (|16p and (|17|l can be 
cast in a matrix form that only deals with the variables d^r 
and dh. Standard relaxation scheme can be applied to the 
resulting matrix. We use uniform grid points in our calcula- 
tions. The grid point number is typically chosen to be 350. 
We implement the radiative boundary conditions such that 
waves propagate away from the density structure in both the 
inner and outer parts of the disk (e.g., Yu & Li 2009). The 
relaxation method requires an initial trial solution that can 
be improved by the Newton— Raphson scheme. After itera- 
tions the initial trial solution converges to the eigenfunction 
of the the two-point boundary eigenvalue problem. 

Figure 2 shows the growth rate of the RWI for differ- 
ent vertical magnetic field strengths of the disc. Real parts 
of the mode frequency for the solid line are all close to 
ujr — 0.99mf2KO. For the dashed line, real parts of the mode 
frequency decrease from ujr — 0.99mS7KO to ujr — 0.95mf2KO 
with the increase of magnetic field. Since Bz is assumed to 
independent of r, the rotation profile of the disc is unaffected 
by Bz when B^ = 0. In this case, we see that as Bz increases. 
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Figure 3. Eigenfunctions of the m = 4 the Rossby modes trapped near the density bump in the disc. The left panels show the non- 
magnetic case, and the right panels have Bz = = 0.78. The top panels show Sur (the solid line is for the real part, the dashed line 
the imaginary part, and the dotted line the absolute value), and the bottom panels show 5u^. 
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Figure 2. The growth rate (in units of Qko) of the m = 4 RWI 
as a function of the plasma /3 parameter [with the corresponding 
vertical magnetic field strength B^ given by Eq. (I20l l] . The solid 
line is the case of B^ = 0, and the dashed line B^ = Bz ■ 



the growth rate first increases slightly and then decreases. At 
(3 ^ 0.3, the growth rate is reduced by more than 10% com- 
pared to the Bz — value. When = Bz, the equilibrium 
rotation profile of the disc is changed from the nonmagnetic 
disc. In this case, the growth rate increases monotonically 
with increasing Bz , reaching about 10% above the nonmag- 
netic value. Note that the dashed hne terminates at /3 ~ 2, 
since an even stronger magnetic field would make nega- 
tive. Comparing the two cases depicted in Fig. 2, it appears 



i' 
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Figure 4. The angular momentum flux associated with the 
Rossby mode. The solid line is for the nonmagnetic case, and 
the dashed line for Bz = B^ = 0.78. In both cases, the eigen- 
function is normalized so that the majcimum Sur equals unity 
(near r = rg). The higher angular momentum flux at r ~ ro for 
the magnetic case is consistent with the faster mode growth rate. 

that the increase in the RWI growth rate mainly results 
from the deeper minimum of the disc equilibrium vortensity 
profile induced by the inclined [B^ > 0) magnetic field. 

Figure [3] gives some examples of the m = 4 eigenfunc- 
tions of Rossby modes trapped around the corotation radius, 
which is also close to the density bump. The left panels show 
the non-magnetized case, and the the right panels show the 
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magnetized case with B+ = = 0.78. Tiie amplitude of 
eigenfunctions are normalized so that the maximum abso- 
lute value of the radial velocity perturbation \5ur\ equals 
unity (this maximum occurs at r ~ ro). The increase of 
5ur near the inner boundary is due to the geometric effects 
(Meheut, Yu & Lai 2012). A Rossby mode excited around 
ro by the density bump can radiate into both the inner and 
outer parts of the disk as spiral density waves. The relative 
phase shift between the real and imaginary parts indicates 
this propagation. 

To understand the origin of the enhanced RWI growth 
rate for magnetized discs, we show in Fig. |4] the angular 
momentum flux associated with the eigenmode for the non- 
magnetic disc and for the disc with — B^ — 0.78. The 
angular momentum flux F(r) across a cylinder of radius r 
is given by (e.g., Goldreich & Tremaine 1979) 



TiSurSu^d(j) y = TrEr^SR (^SurSu'^^ (21) 



F{r) 



where () designates time average and the superscript * de- 
notes complex conjugate. Note that waves carry negative 
(positive) angular momentum inside (outside) the corota- 
tion. The net positive (outward) angular momentum flux 
F{r) around the corotation (close to the density bump) in- 
dicates the growth of the RWI. Higher angular momentum 
fluxes imply higher instability growth rates (see Fig. |4]). 

We have attempted a variational-principle analysis to 
determine the direct effect of magnetic force on the RWI 
(for a given vortensity profile). Such an analysis did not yield 
simple, unambiguous results, and therefore is not presented 
here. This is consistent with the non-monotonic behavior 
exhibited by the solid line in Fig. 2. 



determined by 

o2 IdP B, . B, frrdBA 

To derive the modified perturbation equations including 
5f, it is convenient to define a new perturbation variable 511 
in place of Sh: 

BJB. ^26) 



SP , BJB, 2Sp 

oil = — — ~ c, 

E 47rp ' p 



After some algebra, the final disc perturbation equations can 
be written in the following form: 



dr 



f2mQ 1 
\ rcD r 



rflnE 
dr 



+ 



rr? B, 



(27) 



dr 



511 = [ UJ — K 



B+ 
27rE 



/ 2mQ, 
\ ruj 
2mn B 



where 



rui 27rE 
cl d 



+ D7- 



dlnp 



sn 



ci + cl dr 
1 



In- 



B, 



D7 = 



p2 I p2 ' 

l d{p+^Bl 
p dr 



dr 

2^ dr '• 



" 47rp' 

P 



dr 



(28) 

(29) 
(30) 

(31) 



4 EFFECTS OF FINITE DISC THICKNESS 
4.1 Model Equations 

In Sections 2-3, we have adopted the infinitely thin disc ap- 
proximation [H/r — > 0) and the "internal" disc magnetic 
force is neglected. When finite thickness of the disc is con- 
sidered, an "internal" magnetic force term should be added 
to the right-hand side of Eq. ([2]): 



dz 



(B- Vx)B 



(22) 



Obviously, to include the 3D effect rigorously would require 
examining the vertical stratification of the density and mag- 
netic field inside the disc - this is beyond the scope of this 
paper [see Meheut et al. (2012) and Lin (2012) for the case 
without magnetic field]. Here we consider a simple model 
where the internal density of the disc is assumed to be in- 
dependent of z, so that 



E = 2pH, P = 2pH. 
Then the internal magnetic force simplifies to 



47rp 



B' 



+ (B ■ Vj 



(23) 



(24) 



We assume that only vertical field exists inside the un- 
perturbed disc. The equilibrium rotational profile is then 



4.2 Results 

Finite disc thickness has two effects on the RWI: one is the 
change in the equilibrium rotation profile and the other is 
the direct effect of 5t. We first examine the second effect by 
considering the case of constant Bz (no dependence of r). 
Figure [S] shows the result. We flnd that for the parameter 
considered {ca ~ 0.1, or H/r = 0.1), the RWI growth rate is 
only slightly decreased (by less than 5%) compared to the 
H/r ^0 limit. 

Now consider the effect of the magnetic field gradient. 
With finite disc thickness, such a gradient modifies the equi- 
librium rotation profile. We use the following B^ profile: 



Bz = Bn 



+ -Bamp tanh 



r — ro 
A 



(32) 



Thus, the vertical magnetic field changes from Bn 



at (r - ro)/A < -1 to Bmcan + -Bamp at (r - ro)/A > 
1. The length scale of the Bz variation is the same as the 
size of the density bump. The plasma /3 defined in Eq. H20|l 
is evaluated at r = ro (where Bz = -Bmcan). We assume 
Br+ = Bz- Figure [6] shows some examples of the vortensity 
profiles. The blue curve is for the non-magnetized case, the 
dashed line is for the magnetized case where -Bmoan = 0.64, 
Bamp ~ — 0.6-Bmcan. The magnetic field gradient in Eq. (|25|l 
alters the vortensity profile dramatically. The dashed line 
profile corresponds to the highest growth rate (dotted line, 
where /3 ~ 3) in Fig. [5] 

The dotted line of Figure [5] shows the RWI growth rate 
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Figure 5. Growth rate of the m = 4 RWI as a function of the 
plasma parameter for discs with finite thickness. The solid hnc 
is for the case of constant with = 0, while the dashed 
line for the case of constant B^ = Bz- Note that in these two 
cases, the growth rates are very close to the cases depicted in 
Fig. [2I (where the internal magnetic force 5f is neglected). The 
dotted line shows the case where Bz (r) is given by Eq. Ij32|l with 
Bamp = -0.6-Bmoan and B:^ = Bz- As -Bmoan increases (or /3 
decreases), the growth rate first decreases and then increases. 




Figure 6. Vortensity in both non-magnetized (blue curve) 
and strongly magnetized (red curve, Smcan = 0.64, -Bamp = 
— 0.6-Bmean) situations are shown. The magnetic field gradient 
change the vortensity profile significantly. 



as a function of /3 in the case of Bamp ~ — O.GBmoan- We 
see that the growth rate depends on /3 in a non-monotonic 
way: As /? decreases (-Bmean increases), the growth rate first 
decreases, and then increases, reaching a factor of two com- 
pared to the non-magnetic value (e.g., at /3 ~ 3, correspond- 
ing to the vortensity profile depicted by the dashed line of 
Fig. 6, the mode frequency is LUr = 1.02mf2KO with m = 4, 
and the growth rate is uji — 0.32^x0 )• 

Figure [7] and Figure [8] illustrate how the RWI growth 
rate depends on the Bz gradient, for two different values of 
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Figure 7. The dependence of the growth rate of the m = 4 RWI 
on on the Bz gradient. The Bz profile is given by Eq. I|32| l. with 
Smcan = 0.5 (corresponding to /3 = 5) and B^ = Bz- The solid 
line is for -Bamp < while the dashed line for -Bamp > 0. In such a 
strong magnetic field regime, the magnetic field gradient (regard- 
less of its sign) increases the instability growth rate. Note that 
the lines terminate at the point where the equilibrium rotation 
profile has re'^ = at some radius near the density bump. 
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Figure 8. Same as Fig. [Tlexcept for -Bmoan = 0.1 (corresponding 
to (9 ~ 120). In such a weak magnetic field regime, the growth rate 
of the RWI tends to be reduced by the magnetic field gradient. 

Bmoan- For Small -Bmoan, the growth rate is reduced by the 
Bz gradient (regardless of the sign of he gradient); for larger 
Bmoan, the growth rate is enhanced by the Bz gradient. 



5 DISCUSSION AND CONCLUSION 

In this paper we have carried out linear analysis of the 
Rossby wave instability (RWI) in accretion discs threaded 
by large-scale magnetic fields. The instability can develop 
around the extremum of disc vortensity [(" — K^/(2f2E), 
vorticity divided by surface density], and may play an im- 
portant role in planetesimal formation and angular momen- 
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turn transport in proto-planetary discs, and may also gen- 
erate variabilities in black-hole accretion discs. Our results 
show that the large-scale magnetic field can affect (increase 
or decrease, depending on the field configuration) the RWI 
growth rate even when it has a sub-thermal strength (plasma 
/3 ~ 10 in the disc). For thin discs, the instability growth rate 
can be enhanced by the magnetic field by up to 10%. We 
have also considered discs with finite thickness (but ignoring 
vertical stratofication) , and shown that the RWI growth rate 
can be further increased (by a factor of ~ 2 when /? ~ 1) with 
a steep radial gradient in the magnetic field strength. Just 
like the density bump, such field gradient may be present 
around the transition region between the active and dead 
zones in proto-planetary discs. 

In general, the large-scale magnetic field influences the 
RWI in two ways, either through the direct effect of mag- 
netic force on the perturbed fluid, or by modifying the equi- 
librium disc vortensity profile. We find that the first effect 
is typically smaller than the second. 

Previous studies have shown that an ordered toroidal 
magnetic field inside the disc can suppress the RWI, even 
when the field strength is sub-thermal (Yu & Li 2009). This 
comes about because the in-disc toroidal field changes the 
behavior of wave absorption at the corotation resonance 
(see Fu & Lai 2011). By contrast, the large-scale poloidal 
field considered in this paper does not qualitatively affect 
the corotation resonance, and in most cases leads to an en- 
hanced RWI. Our calculations of discs with finite thickness 
showed that the magnetic field gradient in the radial di- 
rection can significantly affect the dynamics of RWI. More 
rigorous treatments, taking into account the vertical strati- 
fication of the disc and including both poloidal and toroidal 
magnetic field, would be useful. 

The magnetic field - disc configurations studied in 
this paper are susceptible to the production of magneto- 
centrifugal winds. Recent works suggested that in proto- 
planetary discs, MRI-induced turbulence may be inadequate 
for transporting angular momentum and magnetic winds 
may be necessary (Bai & Stone 2011). Our result in this 
paper implies that the RWI can operate under the wind- 
producing conditions, potentially contributing to angular 
momentum transport in the disc. 
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